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SUMMARY 


The application of conventional optimization schemes to aerodynamic design problems 
leads to inner-outer iterative procedures that are very costly. In this report, an alternative 
approach is developed based on the idea of updating the flow variable iterative solutions and 
the design parameter iterative solutions simultaneously. The optimization scheme is suitable 
for application to general aerodynamic problems; here, it is applied to the problem of 
optimizing advanced propeller designs, specifically the SR-3 propeller. The Euler equations 
are assumed to be the flow governing equations in this application, and an implicit 
approximate factorization scheme is used to compute the flow field around the advanced 
high-speed SR-3 propeller. In the computations, the propeller efficiency is maximized subject 
to a given power constraint. The twist distribution of the propeller blade is assumed to be 
described by a polynomial. The coefficients of the polynomial are the design parameters. A 
1.2% increase in the propeller efficiency is achieved. 

Computations were performed to test the scheme’s efficiency, accuracy, and sensitivity. 
The results indicate that the cost of solving an optimization problem with L design parameters 
is approximately equal to L times the cost of solving a regular analysis problem. The scheme 
is highly accurate in determining the solution of the constrained optimization problem. The 
convergence rate of the solution is weakly sensitive to variations in the computational 
parameters and the initial iterative guesses for the design parameters. 


INTRODUCTION 

Solutions of constrained optimization problems minimize an objective function, E, subject 
to given constraints. In aerodynamic applications, the objective function and the constraint 
functions, /,, i = 1, 2, ..., depend on the flow field solution, jf. The optimization scheme 
developed here is applicable to situations in which the flow governing equations are nonlinear 
equations that are solved iteratively. 

Conventional optimization methods (e.g., the steepest descent method and the conjugate 
gradient method) are iterative procedures that require the evaluation of the objective function 
many times before the converged optimum solution is determined. Since E and /, are 
dependent on the flow solution, ~g, in addition to the vector of design parameters, P, the flow 
governing equation must be solved each time E and /,• arc evaluated. Therefore, the 
application of conventional optimization schemes to aerodynamic design problems (refs. 1-5) 
leads to two-cycle (inner-outer) iterative procedures. The inner iterative cycle solves the 
analysis problem for iteratively, while the outer cycle determines the the optimum P 
iteratively. An alternative to this costly procedure is the single-cycle approach, which 



modifies the iterative procedure for solving the flow governing equations so that if and P are 
updated simultaneously (ref. 6). Difficulties have been encountered, however, in attempting to 
apply this approach to advanced propeller design problems (ref. 7). Our objective here is to 
develop a scheme based on the idea of simultaneously updating the flow variables and the 
design parameters that can overcome the problems previously encountered. 

A resurgence of interest in recent years in the turboprop propulsion system has been 
caused by the projected high fuel costs in the 1990’s and the potential savings in fuel con- 
sumption that can be achieved with such a propulsion system. Advanced propellers operate at 
transonic speeds. Therefore, one of the two basic elements required for optimizing the design 
of these propellers is an analysis code capable of solving the nonlinear flow equations about 
the propeller so that the compressibility effects are predicted. The second element is an 
optimization scheme that can be efficiently combined with the analysis code. 

Procedures have been developed for designing propellers by combining vortex lattice aero- 
dynamic analysis methods with standard optimization schemes (refs. 8,9). However, the first 
attempt to optimize propeller designs by using the full potential formulation (ref. 10), which 
includes the necessary elements for transonic design, encountered difficulties in maximizing 
the propeller efficiency subject to a given power constraint (ref. 7). The optimization 
scheme’s inaccurate determination of the constraint surface resulted in these difficulties. 
Thus, in this scheme, efficiency was replaced as an objective function by an approximation, 
which is valid only under special conditions, and computations were limited to low Mach 
numbers. 

In the present work, an optimization procedure is developed based on the idea of updating 
the flow variable iterative solutions and the design parameter iterative solutions simultane- 
ously. This procedure has several common elements with the scheme of Reference 7. How- 
ever, it is more reliable and, thus, eliminates the difficulties encountered by the scheme of 
Reference 7. Although applied here to the propeller design problem, this optimization scheme 
is suitable for application to general aerodynamic design problems. The Euler equations are 
assumed to be the flow governing equations. An implicit approximate factorization scheme 
(ref. 1 1) is used to compute the flow field about an advanced high-speed propeller. 


APPROACH 

The propeller design problem is cast into an optimization formulation in which the 

> * 

optimum design parameter vector, P , is to be determined such that 

E(P ;£*) = min E(P\g) (1) 

p 
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subject to the constraint 


f(P\t) = 0 


( 2 ) 


with the flow variable vector ^ satisfying the flow governing equation 

D(T) = C 


(3) 


subject to the boundary condition 

P) = 0 (4) 

Our objective is to maximize the propeller efficiency, tj. The objective function is therefore 
defined by 

E = - V 

The propeller power requirements are constrained to a specified value through the constraint 
function 


/ — C p *-p o 

Equation (3) is the system of Euler equations governing the flow field, and Equation (4) is the 
propeller solid wall boundary condition. The vector of design parameters, P, defines the pro- 
peller geometrical configuration. 


The goal of the optimization scheme is to determine the values of the design parameters 
that minimize the objective function, E , subject to an equality constraint. A search must 
therefore be conducted in the design parameter space P for the optimum solution, P . This 
optimization problem is most conveniently solved in the rotated design parameter space P, 
with the P x coordinate normal to the constraint surface and the Pj coordinates, where / = 2, 3, 
. . . , L, parallel to the constraint surface. For fixed values of the components of P, let 


_,n+l 

g 


» = o, i, 2, . 


(5) 


be the iterative solution for the analysis problem, where denotes the solution obtained by 
applying the iterative scheme for solving the Euler equations once using g as an initial guess. 
An implicit approximate factorization scheme is used here to solve the Euler equations. It is 
described in Reference 11. As for the analysis solution, obtaining the optimization solution 
requires the repeated application of Equation (5) to update the flow field. While P is held 
fixed in the former case, it is allowed to vary in the latter. The scheme used to update P fol- 
lows. 
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( 6 ) 


The vector of design parameters P is updated every AN iterations. Therefore, 

P n+1 = P n + SP n+1 

where 


SP n+1 = 0, (n+l)/AN* 1,2, 3, ... 

In the iterative steps that satisfy the relation (n + l)/AN= 1, 2, 3, .... the incremental values 
for the design parameters are given by 


where 


*£? +1 = " [min (C m, »>„«)] 


6P? +1 = min 


1, 


SPr 


iA/>r +1 i 


A P|* +1 , l = 2,3,..., L 


f n =/(p n -,g’) 

A P? +1 = - + 1) + c 2 (rr +1 - 1)] 6P? +1 -* n 


(7) 

( 8 ) 


(9) 


A (5P, n+1 - AAr 

T ; +1 

| A £? SPf- +1 ~ AN \ 

A E? = E(? + - E^f*) 

e is a small positive constant and Tf , / = 1, 2, . . . , L, are the set of orthogonal unit vectors 
along the axes of the rotated coordinate system Pf, P$, . . . , P n L . The solution tf? is a solution 
in which the I th component of P is perturbed by e. 

The incremental displacement in the design parameter space introduced so that the con- 
straint may be satisfied is taken in the direction normal to the constraint surface and is deter- 
mined by the chord method in Equation (7). The constant 6P ^ sets an upper limit on the 
magnitude of this incremental displacement. The incremental displacements given by Equa- 
tion (8) are introduced along the coordinate axes, which are parallel to the constraint surface 
with the purpose of reducing the value of the objective function. The sign of the incremental 
correction 6P? +1 , where 6Pf +l is the I th component of the vector 6P n+ \ is chosen to be oppo- 
site to that of dE / dPf. The magnitude of the increment 6P“ +1 is given by 

| <5P( n+1 1 =c | £P ( " +1-AiV | 

with an upper limit given by SP^, where c > 0. If the signs of 6P? +1 and 5P j n+1 - AJV are in 
agreement, then the last two iterative solutions P{* and P ( " _AjV fall to one side of the point along 
the P) direction at which E is a minimum. In this case, c is set equal to the constant c j, which 


4 
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is greater than 1. Increasing the magnitude of the step size in this manner accelerates the 
approach toward the point along the Pj direction at which £ is a minimum. On the other 
hand, if the signs of 6P, n+l and 5P, n+1_AJV are not in agreement, then P" and PP -AJtf fall on oppo- 
site sides of the point along the P, direction at which £ is a minimum. In this case, c is set 
equal to the constant c 2 , which is less than 1. Decreasing the magnitude of the step size in this 
manner is necessary for convergence to the point along the Pj direction at which £ is a 
minimum. 


The updated components of the design parameter vector P* +1 are used to calculate the 

n+1 


new flow iterative solution, g , given by 

g = V’Gf’; 

and the perturbed solutions ifT + \ / = 1, 2, . . . , L, given by 

gl = Mg l, P +€I» ) 


( 10 ) 


( 11 ) 


While the optimization procedure is most suitably conducted in terms of the transformed 
parameters P ( , / = 1, 2, . . . , L, the flow solution is computed in terms of the physical design 
parameters P, , / = 1, 2, . . . , L. To express the transformed design parameters in Equations 
(10) and (1 1) in terms of the original design parameters, it is necessary to use the transforma- 
tion equation, which relates these two sets of parameters. This equation is 

p n+1 = r n+1 p n+1 


where the orthogonal transformation matrix T n+1 is given by 

T'* 1 =[it +1 Tt +l ...7t + ‘l 

The unit vector 7? +1 is normal to the constraint surface at P = P” and is given by 

7? +1 =V/(?"r")/ |V/(P“r")l 


(12) 


where an estimate for C", the I th component of V/, is given by 

ct = [/(P n + «7f; rT) - f(P n ; F“)l / « (13) 


The Gram-Schmidt orthogonalization process, which uses a set of L linearly independent vec- 
tors to construct a set of L orthonormal vectors, is used to construct the unit vectors 

7f +1 , / = 2, 3 L, along the rotated axes P, n+1 , / = 2, 3 L. The following equation is 

used for this purpose: 


-/»+! 

h 


-^•n+1 

ITT 1 ! 


/ == 2, 3, . . . 
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where 


v +i = Tf-Eef-Tr 1 )^ 1 

k=i 

In the initial iterative step, the vectors T, are given by if =T,,/= 1, 2, .... L, where 

1 1 1*2 , . . . , L, are the set of orthogonal unit vectors along the axes of the coordinate sys- 

tem P u P 2 , P L , 

While the flow variable vector f is updated each iterative step, the coordinate system in the 
design parameter space is rotated every AN iterations. The unit vectors 7J, like the vector of 
design parameters P, are updated only in the iterative steps that satisfy the relation 
(n + l)/AN = 1, 2, 3, ... . 

The optimization scheme described above requires that L + 1 iterative problems be solved 
in parallel. In addition to the main solution, L perturbed solutions are computed in which 
each of the design parameters in the transformed space P u P 2 , . . . , P L is perturbed. The 
computational costs and the computer memory requirements are therefore proportional to 
Z. + 1. A modification to this scheme, which requires that only L iterative solutions be 
obtained, is now introduced. In the modified procedure, the perturbation solution associated 
with the perturbed design parameter in the direction of the Pi axis, normal to the constraint 
surface, is not computed. This solution was used in Equation (13) to compute GJ, which is 
required for the calculation of the vector 7? +1 , which determines the direction normal to the 
constraint surface in Equation (12). In the absence of this solution, a new procedure for rotat- 
ing the design parameter space must be defined. The procedure is first explained for the case 
of a two-design-parameter problem, and then it is extended to the general multi-design- 
parameter problem. 


Figure 1 shows the design parameter space for a two-design-parameter problem. In the fig- 
ure, the constraint function values /J /?, /£ are defined as follows: 


f n o = f(P\t n ) 

n = f{T + ti) 
/ 2 » = /(F + £ t£^) 


In the modified procedure, the chord method, used in Equation (7) to satisfy the constraint 
condition, is used to rotate the design parameter space. The rotation angle SOty 1 given by 


6dtf l = tan" 1 


c(/ 2 n - /;) 

6 


(14) 


is used to rotate the coordinate system, where the subscript M indicates that the modified 
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pn+1 
— 2M 


\ 

\ 

pn+ 1 



Figure L Two-Dimensional Design Parameter Space 


scheme is used. The angle is now compared to the corresponding rotation angle 60 n + l 
used in the original scheme and given by 

66 n+ 1 = tan' 1 h—Lz. i ( 15 ) 

« n-n\ 

This comparison shows that the term /“ — /J in the original scheme is replaced by t/C in the 
modified scheme. Therefore, the modified scheme may be viewed as the original scheme with 
the exception that the exact value forTi is replaced by an approximate estimate in which the 
gradient of / in the direction of the axis, Cj, is not calculated but is estimated using the 
same proportionality constant used in the chord method of Equation (7). Thus, 

= £ (16) 

This is applicable for both the two-design-parameter problem and the general multi-design- 
parameter problem. 
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In the optimization scheme developed here, corrective increments are applied to the design 
parameter solutions every few iterations of updating the flow solutions. For convergence to 
occur, the signs of the increments must be chosen correctly to allow the iterative solution to 
approach the desired solution. The magnitudes of the increments are dependent on the com- 
putational constants c j, c 2 , and C. Because the design parameters are updated frequently dur- 
ing the iterative process, we are not concerned with determining the incremental step sizes 
that lead to the highest short-term convergence rate. In fact, this may be difficult to define, 
since the flow variable solutions are continuously changing during the iterative process. Our 
aim is to achieve design parameter convergence over a long term defined by the number of 
iterations required for the flow solution convergence. A wide range of incremental step sizes 
should produce the desired convergence properties over many iterations, even though conver- 
gence properties over a few iterations may differ. These comments apply to both of the 
schemes described above for determining the design parameter space rotation. The direct 
procedure for determining the design parameter space rotation in the original scheme is 
replaced by an iterative procedure in the modified scheme. Since this rotation is updated fre- 
quently during the iterative process, this replacement should have no substantial effect on the 
overall convergence of the solution. 

A potential problem exists when the modified scheme is used for rotating the design 
parameter axes. This problem is now discussed, then suggestions for overcoming it are 
presented. 

In the first AN—\ iterative steps of solving the problem, the coordinate system in the 
design parameter space coincides with the original unrotated design parameter space 
P 1 , P 2 , ■■■ , Pl ■ At the AN 01 iterative step, a new rotated coordinate system is determined. 
When Equation (13) for determining Cf* -1 is used, we are guaranteed that the vector lf N 
points in the direction in which the constraint function increases. Consequently, the use of 
Equation (7) will cause the iterative solution to approach the constraint surface. When Equa- 
tion (13) is replaced by Equation (16) for determining there is a possibility that the 

computed vector 7^ will point in the direction in which the constraint function decreases. In 
this case, the assumption that C is positive is wrong, and using it will cause the solution to 
diverge. This occurs if the vector '? 1 is nearly in the direction of — V/^ -1 . That is, if the 
quantity 

y yAJV-i 

“iV/AW" 1 ) *1 

is close to unity. The probability of this occurring is approximately 1:4 in a two-design- 
parameter problem and is reduced further as the number of design parameters increases. 
There are two suggested approaches for overcoming this problem. In the first approach, the 
initial few iterations are performed using the original scheme for determining Cf by Equation 
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(13) in order to determine the correct initial directions for the Pj axis. This may then be 
updated using the modified scheme. Equation (16) in the rest of the computation. Realizing 
that the probability for the potential problem to occur is small, the second approach uses the 
modified scheme from the beginning of the computation. If divergence does occur, then the 
constraint function is redefined to be equal to the negative of the original constraint function, 
and the problem is solved again. 


RESULTS 


The optimization procedure described above, combined with the Euler analysis code 
developed by Yamamoto et al. (ref. 11), was used to find the twist distribution for the blades 
of the eight-bladed SR-3 propeller with the objective of maximizing its efficiency under the 
constraint of a desired power coefficient given by C,*, = 1.7. The computations were per- 
formed for a free-stream Mach number of 0.8 and an advance ratio of 3.06. We let 0 o »/ t be the 
blade angle at the 75% blade span corresponding to the desired power coefficient, and we 
took the blade angle distribution, p o (r), corresponding to this propeller as our base configura- 
tion. A perturbation, y9'(r), to the blade twist distribution, 0 o (r), was computed so that the 
propeller efficiency is maximized subject to the power constraint. The perturbation twist dis- 
tribution is given by 


fi\r) = Pi + P 2 


r-R/2 
R/ 2 


+ 2 P 3 


r-P/2 

2 f 

P/2 

2 


(17) 


where Pi, P 2 , and P 3 are the components of the vector of design parameters P and R is the 
propeller radius. 


Experimentation with the propeller analysis code indicated that the flow iterative solution 
diverges when the blade tip angle exceeds a certain limit. To exclude the region leading to 
the divergence from our search in the design parameter space, the following redefinition of 
the objective function was introduced: 

E = -n + max [ 0.0, 0. 1 [v^a+Pl-p) ] 08) 

where /x determines the allowable search region. As the value of /x increases, the allowable 
search region also increases. The value of /x was taken to be equal to 5.0 unless otherwise 
specified. 

The mesh used in the following computations consists of 45 points in the axial direction, 
21 points in the radial direction, and 11 points between adjacent blades in the circumferential 
direction. Computations are initialized by the SR-3 flow solution, which corresponds to a 
54.9* angle at the 75% blade span. This initial solution was intentionally chosen not to be a 
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close approximation of the desired solution. In all the following computations, the modified 
coordinate rotation scheme, which determines Ci by Equation (16) instead of Equation (13), 
is used unless otherwise specified. Also, unless otherwise specified, the initial iterative 
guesses for the design parameters are set equal to zero and the computational parameters Ci, 
c 2 , C, SPl 6P%, 6P £, and AN are given, respectively, by 1.2, 0.6, 3.0, 0.5, 0.5, 1.0, 0.0001, 
and 40. The computations were performed on the NASA Lewis Cray X-MP computer. 

The optimization scheme was mainly tested using the analysis code developed by 
Yamamoto et al. (ref. 1 1) in its original form, referred to here as code O. Towards the end of 
the present study, an error was discovered in the portion of the analysis code that computes 
the propeller performance. This error was corrected, and the resulting modified code is 
referred to here as code M. Some optimization results were recomputed using code M. Both 
sets of computations resulted in substantially different solutions. The convergence properties 
associated with both codes also showed substantially different behavior. While there may be 
no interest in the first set of solutions for the purpose of improving the propeller design, both 
the first and second sets of computations are of equal interest for the purpose of testing the 
optimization scheme. Thus, the results of both sets of computations are presented below. The 
results obtained by using code O are presented first, followed by those obtained by using code 
M. 


Using code O, it was determined that p o > /t = 58.067*. The value of C p for the initial flow 
solution, which corresponds to a /9y„ value of 54.9*, was 1.1. The optimization procedure was 
applied to two-design-parametcr problems and to three-design-parameter problems. For the 
two-design-parameter computations, the values of P 3 in Equations (17) and (18) are set equal 
to zero. Results for the two-design-parameter problem are presented, followed by those for 
the three-design-parameter problem. 

The design parameters predicted by the optimization scheme are given by P[ = -2.83°, 
P 2 = 5.51*. The predicted solution does satisfy the power constraint. The value of C p 
corresponding to this solution is 1.6999. The objective function, E, was reduced from the 
value -0.839 in the case of the original design, with = P 2 = 0.0, to the value -0.908 in the 
case of the optimized design. The value of the efficiency was increased from 0.839 for the ori- 
ginal design to 0.9 1 0 for the optimized design. 

The iterative histories of the design parameters are shown in Figure 2, while the iterative 
histories of the power and efficiency are shown in Figure 3. From these figures two distinct 
stages in the convergence process of the solution may be identified. In the first stage, rela- 
tively rapid changes in the values of P, C p and r) occur as they approach the converged values 
of the solutions. At the end of this stage, these parameters are close to their final values. In 
the second stage, minor adjustments take place as the parameter solutions converge to their 
final values. 
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NUMBER OF (Ti RATIONS 

Figure 2. Design Parameter Iterative Histories for the Two-Design-Parameter 
Problem (Code O) 



NUMBER OF ITERATIONS 

Figure 3. Power and Efficiency Iterative Histories for the Two-Design- 
Parameter Problem (Code O) 



NUMBER OF ITERATIONS 

Figure 4. Residual Iterative Histories for the Analysis Problem and the 
Two-Design-Parameter Optimization Problem (Code O) 
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The residual, R& is a measure of the convergence of the flow field solution. Figure 4 com- 
pares the residual history for the design problem, in which ^ is updated in addition to P, to 
the residual history for the regular analysis problem, in which only is updated while P is 
held fixed. The figure indicates that modifying the propeller geometry in the design problem 
as the iterative solutions for the flow variables are updated does not negatively affect the rate 
of convergence of the flow field solution in comparison to the analysis problem. In fact, the 
following results of our computations show that the convergence of the flow field solution is 
accelerated when the design parameters are updated to satisfy the power constraint or to 
satisfy the the conditions of the optimization problem. For a regular analysis problem with P 
set equal to P , the number of iterations required for convergence was 4710. Throughout this 
report, convergence is assumed to be achieved when the magnitude of the residual, Re, is 
reduced to the value of 10 7 . For a constrained solution in which the second component of 
the design parameter vector, P 2 , was set equal to the value P 2 , while the first component was 
updated throughout the iterative process so that the constraint C p = C ^ would be satisfied, 
convergence was attained after 4040 iterative steps, indicating an increased convergence rate 
relative to the regular analysis problem. For the design problem in which both P x and P 2 
were updated in a manner that allows the constraint C p = C ^ to be satisfied and the objective 
function E to be minimized, the number of iterative steps required for convergence was 
further reduced to 3250. 

On the average, 0.972 cpu second was required for the iterative step in the design problem, 
while 0.403 cpu second was required for the iterative step in the analysis problem. Therefore, 
the average design iterative step required slightly more than double the cpu requirements for 
the analysis iterative step. In the design problem, two analysis problems are solved in parallel. 
The additional cpu requirement for the design problem is mainly due to generating a new 
computational mesh whenever the design parameters are updated. 

For a regular analysis problem, the computational mesh is generated only one time at the 
beginning of the computation. For a design problem, however, it is necessary to regenerate 
the computational mesh whenever the design parameters are updated. In the present compu- 
tations, this was done once every 40 iterative steps. The cost of mesh generation relative to 
the cost of solving the flow equations was acceptably low. As the value of AN decreases, 
however, a point may be reached at which the cost of generating the mesh becomes exces- 
sively high, and it may represent a substantial fraction of the total computational cost. In this 
case, a possible alternative to regenerating new meshes, whenever the design parameters are 
updated, is the use of approximate meshes that are generated by linearly combining L + 1 
reference meshes. The reference meshes may be updated every K AN iterative steps, where 
^ > !• The need for making this approximation does not arise here, as the propeller analysis 
code used here has relatively slow convergence properties and, therefore, the appropriate AN 
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value is relatively large. However, the use of accelerating schemes, such as the multigrid 
scheme, would allow the AN value to be sufficiently low to require the use of the mesh 
approximation discussed above. 

We have performed a single computation using the exact formulation for calculating Ti, as 
given by Equation (12), with computed by Equation (13). This formulation requires solving 
L+l problems in parallel instead of L problems, in the case of the approximate formulation 
given by Equation (16). The average iterative step for this computation required 1.474 epu 
second. The number of iterations required for convergence was 3425. Comparing these values 
to the corresponding values for the approximate formulation indicates that there is a strong 
advantage in using the approximate formulation over the exact formulation. 

To verify that the computed solution is indeed the optimum solution, solutions were com- 
puted that were slightly perturbed from the optimum predicted solutions but that satisfied the 
power constraint. Table 1 compares the values of the objective function for the solution 
predicted by the optimization scheme, shown in the: first row, to those for the perturbed solu- 
tions, shown in the second and third rows. It is apparent from the table that perturbing the 
design parameters causes the value of the objective function to increase. Therefore, the 
design parameters predicted by the optimization scheme do indeed minimize the value of the 
objective function. 


Table 1. Objective Function at Optimum Solution and 
Perturbed Solutions for Two-Design-Parameter Problem (Code O) 


Pi 

P 2 

E 

-2.83 

5.51 

-0.90773 

-2.73 

5.31 

-0.90730 

-2.93 

5.71 

-0.90728 


The sensitivity of the scheme’s convergence to the initial iterative guesses of the solution 
and to the computational parameters was tested by recomputing the problem defined above 
with perturbed initial conditions and computational parameters. Table 2 shows the number of 
iterative steps, n c , required for convergence when different values are used for the initial 
iterative solutions and the computational parameters. It is clear from the table that the con- 
vergence properties of the scheme are weakly sensitive to the values of the initial conditions 
and the computational parameters. Needless to say, there is an optimum set of values for these 
parameters that maximizes the convergence rate of the scheme for a given problem. However, 
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Table 2. Effect of Perturbing Initial Conditions and Computational Parameters 
on Scheme’s Convergence for Two-Design-Parameter Problem (Code O) 


Pi 

r% 

An 

C 

Cl 

c 2 

n c 

0.0 

0.0 

40 

3.0 

1.2 

0.6 

3250 

3.0 

-5.0 

40 

3.0 

1.2 

0.6 

3690 

0.0 

0.0 

25 

3.0 

1.2 

0.6 

3376 

0.0 

0.0 

40 

4.5 

1.2 

0.6 

3252 

0.0 

0.0 

40 

6.0 

1.2 

0.6 

3250 

0.0 

0.0 

40 

3.0 

1.5 

0.6 

3333 

0.0 

0.0 

40 

3.0 

1.2 

0.4 

3120 

0.0 

0.0 

40 

3.0 

1.5 

0.4 

3281 


within a relatively wide range of these parameter values, good convergence is achieved. This 
is due to the frequent updating of the design parameters in the course of solving the problem. 
The cpu requirement for the average iterative step is approximately the same for all the cases 
solved, except for the case in which AN = 25. The cpu requirement for the average iterative 
step in this case is given by 1.078 seconds, in comparison to approximately 0.972 second for 
the other cases. This is due to the increased frequency of generating the computational mesh 
in the case with AN = 25. Figures 5 through 7 show the iterative histories for P u P 2 , rj, C p 
and R e for the case in which the initial iterative guesses for the design parameters, P? and P%, 
were perturbed. Among all the perturbed computations, the rate of convergence for this case 
was affected the most. 

The computations performed above for the two-design-parameter problem were performed 
with a value of 5.0 for /i. To perform computations that allow both parabolic and linear 
modifications to the blade angle distributions, it was necessary to reduce the value of n to 4.0. 
The three-design-parameter optimization computations were solved using this value for \i. 
The main two-design-parameter computation was also repeated using this value for n to allow 
a comparison between the two-design-parameter and the three-design-parameter results. 

The optimum values of the design parameters for the two-design-parameter problem with 
H = 4.0 were found to be given by P[ = -2.35" and P 2 = 4.56". The value of C p corresponding 
to this solution is 1.6999. The objective function E was reduced from the value —0.839 in the 
case of the original design, with P 1 = P 2 = 0.0, to the value -0.897 in the case of the optim- 
ized design. The value of n was increased from 0.839 for the original design to 0.900 for the 
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Figure 5. Design Parameter Iterative Histories for the Two-Design-Parameter 
Problem with Perturbed Initial Conditions (Code O) 



Figure 6. Power and Efficiency Iterative Histories for the Two-Design- 
Parameter Problem with Perturbed Initial Conditions (Code O) 



Figure 7. Residual Iterative Histories for the Analysis Problem and the 

Two-Design-Parameter Optimization Problem with Perturbed Initial 
Conditions (Code O) 
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optimized design. As expected, the magnitudes of both E and 17 determined with M = 4.0 are 
less than those determined with fi = 5.0. As the value of n decreases, the restriction on the 
allowable search region in the design parameter space increases. In the two-design-parameter 
problem, 3235 iterative steps were required for convergence. The cpu requirement per itera- 
tive step was 0.972 second. The optimum values of the design parameters for the three- 
design-parameter problem with /i = 4.0 were found to be P[ = -2.77*, P 2 = 4.50°, and 
P 3 = -1.20°. The corresponding values of C ^ E, and ij are given by 1.6999, -0.900, and 
0.905, respectively, indicating a superior design to that achieved by using only two design 
parameters. The number of iterative steps required for convergence was 3228, while the cpu 
requirement per iterative step was 1.459 seconds. The iterative histories for P u P 2 , P 3 , C p 
and R e are shown in Figures 8 through 10. 

To verify the accuracy of the computed solution, several solutions were computed that 
were slightly perturbed from the optimum predicted solution but that satisfied the power con- 
straint. Table 3 compares the values of the objective function for the solution predicted by 
the optimization scheme, shown in the first row, to those for the perturbed solutions shown in 
the following rows. It is apparent from the table that perturbing the design parameters causes 
the value of the objective function to increase. Therefore, the design parameters predicted by 
the optimization scheme do indeed minimize the value of the objective function. 


Table 3. Objective Function at Optimum Solution and 
Perturbed Solutions for Three-Design-Parameter Problem (Code O) 


Pi 

P 2 

Pz 

E 

-2.77 

4.50 

-1.20 

-0.90026 

-2.87 

4.50 

-1.45 

-0.90011 

-2.87 

4.69 

-1.20 

-0.89986 

-2.67 

4.50 

-0.93 

-0.90012 

-2.67 

4.30 

-1.20 

-0.89988 


The initial attempt to determine the optimum solution for the two-design-parameter prob- 
lem by using code M led to an unconverged solution. An oscillatory iterative history was 
observed for the design parameters (see Figure 11). This differs substantially from the 
corresponding iterative history associated with the use of code O, as seen by comparing Fig- 
ures 2 and 11. The iterative histories for rj, C p and R E obtained in this initial attempt are 
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Figure 8. Design Parameter Iterative Histories for the Three-Design-Parameter 
Problem (Code O) 



Figure 9. Power and Efficiency Iterative Histories for the Three-Design- 
Parameter Problem (Code O) 


0.0 



Figure 10. Residual Iterative Histories for the Analysis Problem and the 
Three-Design-Parameter Optimization Problem (Code O) 
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Figure 11. Design Parameter Iterative Histories for the Two-Design-Parameter 
Problem with = 1.2 (Code M) 



Figure 12. Power and Efficiency Iterative Histories for the Two-Design- 
Parameter Problem with c i = 1.2 (Code M) 



Figure 13. Residual Iterative History for the Two-Design-Parameter 
Optimization Problem with Ci « 1.2 (Code M) 
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shown in Figures 12 and 13. To isolate the effects of satisfying the constraint condition from 
the effects of minimizing the objective function on the observed behavior, a series of tests was 
conducted. The right-hand side of Equation (8) was set equal to zero in some of the tests, 
while the right-hand side of Equation (7) was set equal to zero in others. These tests indicated 
that the convergence properties of codes O and M were similar when the design parameters 
were modified with the objective of satisfying the constraint only; however, they differed 
when the design parameters were modified with the objective of minimizing the objective 
function only. In this latter case, the tests showed that the efficiency computed by code M 
responds to changes in the design parameters at a much slower iterative rate than that associ- 
ated with code O. 

In the computations performed by using code O, the incremental step sizes of the design 
parameters were reduced by the factor frequently enough to allow convergence to occur. 
On the other hand, when code M was used, the frequency at which the incremental step sizes 
were reduced by the same factor was not sufficiently high to allow convergence to occur. A 
satisfactory solution to this problem was obtained by requiring the design parameter incre- 
mental step sizes to decrease continuously as the number of iterations increases. The function 
of the factor c u which was originally used as an incrementing factor with a value greater than 
unity, was therefore changed to a decrementing factor with a value slightly less than unity. 
The value of c x used in computing the following results was 0.98. 

Using code M, it was determined that /9 Q y 4 = 57.648*. The value of C p for the initial flow 
solution, which corresponds to a p>/ t value of 54.9°, was 1.2. The results of three computations 
with different initial iterative guesses for the design parameters are shown in Table 4. 


Table 4. Comparison of Computational Results 
for Different Initial Iteratbe Guesses (Code M) 


/>? 

BIB 

mm 

P 2 

■a 


n c 

0.0 

0.0 

-0.67 

1.22 

1.6998 

0.82336 

3260 

-2.0 

-2.0 

-0.70 

1.27 

5.6998 

0.82331 

4000 

2.0 

2.0 

-0.60 

1.10 

1.6999 

0.82333 

4300 


By optimizing the blade shape, the value of the efficiency was increased from 0.8229 for 
the original design to 0.8233 for the optimized design. The table shows that the computed 
solutions satisfy the power constraint to a high degree of accuracy. The values of r\ predicted 
by the different computations are in agreement to within 5x10“®. Differences among the 
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values computed for the optimum solutions Pi and P 2 are small; however, they are larger than 
those associated with C p and »?. This indicates that the sensitivity of the propeller perfor- 
mance to variations in the design parameters is relatively weak. Further reductions in the 
differences among the results of the three computations may be obtained by continuing the 
iterative computations further. For a regular analysis problem with P set equal to P\ the 
number of iterations required for convergence was 4320. A comparison of this number with 
the number of iterations required to solve the optimization problem, given in the table, shows 
that the cost of solving the optimization problem is approximately twice the cost of solving a 
regular analysis problem. The iterative histories for P u P 2 , n, C p and R E for the computation 
with initial iterative conditions given by F? = F§ = 0.0 are shown in Figures 14 through 16. 

Figure 1 7 shows that the constraint (constant power) curve and the constant efficiency 
contours in the design parameter space are nearly parallel. Along the constraint curve, the 
variation in the objective function is very small, while its variation is relatively large in the 
direction normal to the constraint curve. This configuration causes difficulties when conven- 
tional optimization schemes are used. The success of the present scheme in solving this prob- 
lem is an indication of its reliability. 

In the computations presented above, the effect of varying the linear term of Equation (17) 
on the propeller efficiency was investigated. To investigate the effect of varying the quadratic 
term in Equation (17) on the propeller efficiency, a computation was performed in which F 3 
was allowed to vary while F 2 was set equal to zero. In this case, the design parameters 
predicted by the optimization scheme were given by P[ = -0.19°, F 3 = -2.07°. The value of 
C p corresponding to this solution was 1.7000, and the value of 17 was 0.82549. The number of 
iterations required for convergence was 3980. A comparison of the values of tj for the two 
cases in which (P u P 2 ) and (Pi, F 3 ‘) were the design parameters shows that the introduction 
of a quadratic perturbation to the twist distribution is more effective in increasing the effi- 
ciency than the introduction of a linear perturbation. 

Finally, the optimum values of the design parameters for the three-design-parameter prob- 
lem were found to be Pi = -3.34°, F 2 = 3.92°, and P 3 = -3.23°. The corresponding values of 
C p and ij are given by 1.7000 and 0.83291, respectively. It is apparent that using a combina- 
tion of linear and quadratic perturbations in the blade angle distribution is much more effec- 
tive for improving the efficiency than using only one of these distributions. Relative to the 
original SR-3 design, using both perturbed distributions increased the propeller efficiency by 
0.0100. This is compared to a value of 0.0026 for the quadratic distribution alone and a value 
of 0.0004 for the linear distribution alone. The number of iterative steps required for conver- 
gence was 4380 in comparison to 4460 for the regular analysis problem. The iterative his- 
tories for P u P 2 , F 3 , n, C p and R E are shown in Figures 18 through 20. 
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Figure 14. Design Parameter Iterative Histories for the Two-Design-Parameter 
Problem with C\ = 0.98. (Code M) 



Figure 15. Power and Efficiency Iterative Histories for the Two-Design- 
Parameter Problem with c x = 0.98 (Code M) 



Figure 16. Residual Iterative Histories for the Analysis Problem and the 

Two-Design-Parameter Optimization Problem with c \ - 0.98 (Code M) 
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Figure 18. Design Parameter Iterative Histories for the Three-Design- 
Parameter Problem with C\~ 0.98 (Code M) 



NUMBER OF mi RATIONS 

Figure 19. Power and Efficiency Iterative Histories for the Three-Design- 
Parameter Problem with Ci * 0.98 (Code M) 



Figure 20. Residual Iterative Histories for the Analysis Problem and the 

Three-Design-Parameter Optimization Problem with c x = 0.98 (Code M) 
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Figure 21 compares the optimum blade angle perturbations from the SR-3 baseline design 
predicted for the cases of linear, quadratic, and combined linear and quadratic shape func- 
tions. Curve C, which gives the blade angle perturbation distribution for maximum improve- 
ment in efficiency, shows that the efficiency of the SR-3 propeller can be improved by reduc- 
ing the blade angle distribution both at the hub and at the tip. This explains the observed 
weak sensitivity of the propeller efficiency to linear variations in the blade angle distribution. 
The use of a linear shape function allows an increase in the blade angle at either the tip or the 
hub positions and a decrease in the blade angle at the other position. Therefore, the positive 
effect on efficiency resulting from the perturbed blade angle distribution at one of these posi- 
tions tends to cancel the negative effect resulting from the perturbed blade angle distribution 
at the other position leading to the apparent insensitivity of the efficiency to linear variations 
in the blade angle distribution. The maximum improvement in efficiency obtained here 
resulted from the use of linear and quadratic shape functions. Further improvement may be 
obtained by using other shape functions. 



Figure 21. Optimum Blade Angle Perturbations. (Code M) 
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CONCLUSIONS 


In this report, we developed a scheme for solving constrained optimization problems in 
which the objective function and the constraint function are dependent on the solution of the 
nonlinear flow equations. The scheme updates the design parameter iterative solutions and 
the flow variable iterative solutions simultaneously, thereby eliminating the need for the costly 
inner-outer iterative procedure associated with the use of conventional optimization schemes. 
The scheme was applied to the problem of optimizing an SR-3 advanced propeller design with 
the Euler equations assumed to be the flow governing equations. The optimized design caused 
a 1.2% increase in the propeller efficiency. 

Computations were performed to test the scheme’s efficiency, accuracy, and sensitivity. 
Two computer codes were used in the computations. The results from both codes indicate 
that the cost of solving an optimization problem with L design parameters is approximately 
equal to L times the cost of solving a regular analysis problem. The scheme is highly accurate 
in determining the solution of the constrained optimization problem. The scheme’s sensitivity 
to the computational parameters was tested using only one of the two codes. This study 
showed that the convergence rate of the solution is weakly sensitive to variations in the com- 
putational parameters and the initial iterative guesses for the design parameters. 
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APPENDIX: NOMENCLATURE 


Cl 

Cl 

C 

C P 

r 

v -'pO 

D 

e't 

E 

f 

f 

G, 

U 

it 

L 

«c 

P 

P 

Pi 

Pi 

r 

R 

Pe 

P'U 

PoH i 

Po 

P‘ 

sP 

^P max 

AN 

€ 

V 

p 

p> 


incrementing factor for optimization scheme [see Equation (9)] 

decrementing factor for optimization scheme [see Equation (9)] 

positive constant for chord method [see Equation (7)] 

power coefficient 

desired power coefficient 

propeller diameter 

unit vector along the P[ axis 

objective function 

constraint function 

solution of the flow governing equations 

I th component of V/ relative to rotated c oordinate system 

unit vector along the Pi axis with components defined relative to the unrotated 

design parameter coordinate system 

unit vector along the Pi axis with components defined relative to the rotated 

design parameter coordinate system 

number of design parameters 

number of iterations required for convergence 

vector of design parameters 

vector of design parameters relative to rotated coordinate system 
I th component of design parameter vector 

I th component of design parameter vector relative to rotated coordinate system 

radial coordinate 

blade tip radius 

residual Euclidean norm 

SR-3 blade angle at 75% blade span 

P> /t which corresponds to the power coefficient C ^ 

unperturbed blade angle distribution 

blade angle distribution perturbation 

incremental vector used to update the vector of design parameters 
maximum incremental value allowed in updating the design parameters 
number of iterative steps at which P is periodically updated 
small positive incremental value used to perturb the design parameters 
efficiency 

parameter determining the allowable region in design parameter space for 
searching for the optimum solution [see Equation (18)] 
flow iterative solution 
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APPENDIX: NOMENCLATURE (CONT.) 


Superscripts 

n iteration number 

* optimum value 


Subscripts 

M coordinate system rotated by the modified scheme 

_ rotated coordinate system 
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